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Abstract 
Introduction. Long-stroke movements in automated pneumatic drives account for a significant number of executive 


movements in coordinate tables, automated warehouses, cutting machines, etc. Long-stroke movements degrade the 
dynamic quality and positioning of the drive. This is due to the friction of the piston and the nonlinear characteristics of 
the compressed gas flow in significant volumes of the pressure and drain cavities of the cylinder. Thus, it seems 
promising to create an automated position pneumatic actuator for long-stroke movements. This will increase the 
productivity of processes while providing the declared accuracy. The objective of the work is to obtain a mathematical 
model and dependences of the critical parameters of the proposed automated position long-stroke pneumatic drive of 
fabrication system in the areas of acceleration, steady-speed movement, deceleration, and braking. 

Materials and Methods. The basis for calculations and modeling was the scheme of two trajectories of movement from 
point A to point E, taking into account the forces expended on these processes. The optimal displacement was 
determined using the Portnyagin’s principle (i.e., optimal performance). Proportional drive control was presented as a 
method of achieving the result. For long-stroke drive movements, schematic solution and design scheme were 
visualized in detail (presented as drawings). An original jet sensor with an internal pneumatic connection and a 
pneumo-mechanic discrete-proportional device for the control loop performance were proposed. The mathematical 
model included the movement and braking of the piston, the balance of mass flow, the pressure at points, and the 
control loop. The system of equations was solved by the Runge-Kutta method in the SimInTech software product. 
Based on the results of the study of a generalized mathematical model, the dependences of changes in the kinematic, 
power and pneumatic properties of the drive were constructed in real time during a typical positioning cycle. The 
information was summarized and presented as a set of graphs. 

Results. The mathematical model was formed according to a set of calculations. It took into account the dependences 
characteristic of the movement of the piston of the pneumatic cylinder. The balance of mass flow was investigated by 
the equations of gas flow during compression in the chamber, through distributors and chokes, in the discharge and 
drain cavities and in the control device. Inequalities describing the pressures at the points and the control loop were 
considered. A complex mathematical model was solved in the SimInTech software environment by the Runge-Kutta 
method with a variable integration step. A fragment of the program was selected as one of the illustrations. It showed 
that the software used the following indicators for calculations: target and reduced coordinates; absolute gas constant; 
coefficients of spring stiffness, resistance, adiabatic and viscous friction in the piston; compressor pressure; mass of the 
moving parts of the pneumatic actuator; strength of external resistances; diameters of the pipeline, the pneumatic 
cylinder piston and the braking device; length of the stroke of the cylinder piston; area of piston cavities and throttles; 
length of the pipeline and its internal volume. Thus, the program manipulated a significant set of data, which made it 


possible to obtain meaningful and adequate results. The relationship of blocks and diagrams used in solving the model 
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was schematically shown. We are talking about graphs of movements, areas, pressures, velocities and temperatures. 
Blocks with the program text and intended for integration were used. Thus, a mathematical model of an automated 
pneumatic drive of the fabrication system and the dependences of the basic parameters of its operation were obtained. 
The graphs indicated that the operating mechanism of the pneumatic actuator properly followed the proposed trajectory. 
Discussion and Conclusion. The research results allowed us to consider several stages of long-stroke movement of the 
drive, to determine the time frame of these processes (from 0 to 0.65s), as well as changes in pressure and speed of 
movement of the pneumatic cylinder carriage recorded in these intervals. There were five such stages: acceleration, 
steady-speed movement, deceleration, movement with positioning speed, and braking. Further research will focus on 


optimizing the system to reduce the duration and maintain accurate positioning under external influences. 
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AHHOTalna 

Beedenue. Ha mHHOXOOBEIe NepeMeLIeCHHA B ABTOMATH3MPOBaHHbIX THEBMONPHBOaX MPHXOAMTCA 3HAYMTEIIbHOe 
KOJIM4eECTBO HCHOJIHUTCIbHBIX JBHKCHHM B KOOPAMHAaTHbIX CTOIAaX, Ha ABTOMATH3HPOBaHHbIX CKIayax, pacCKpOHMHBIX 
MalmMHax UT. 1. J MHHOXOJOBbIC NepeMeleHHA YXyMMaloT WHHAaMMYeCKOe KayeCTBO H NO3HIMOHUpOBaHHe MpuBora. 
3ro obycioBieHo TpeHHeM MOpWHA HU HeMHeMHbIMM XapaKTepHCTHKaMH TOTOKa CxKaTOrO ra3a B 3HAYMTeJIbHBIX 
oObeMax HanopHoli U CIMBHOM MoNocTeHi WuMH Apa. Takum oOpa30M, ipecTaBsiaeTCA MepCHeKTHBHBIM CO3jaHve 
AaBTOMATH3HPOBaHHOTO NO3HIMOHHOTO MHEBMOMPHBOAa WIA JVIMHHOXOAOBBbIX TepeMelleHHH. ITO MO3BOJIMT TOBbICHTb 
TIPOW3BOMTeIbHOCTb TpoyeccoB mpu obOecneyeHHH 3aaBIeHHOH TouHocTH. Ilemb padoTbl — nomy4enne 
MaTeMaTHYeCcKOH MOeIM U 3aBHCHMOCTeli OCHOBHBIX THapaMeTpOB [peIO%KeHHOTO aBTOMATH3HpOBaHHOrO 
MO3HUMOHHOFO JVIMHHOXOAOBOrO MHEBMOMpHBOa TeEXHOOrMYecKorO OOOpyAOBaHHA Ha y4YacTKax pa3roHa, IBWKCHHA 
C YCTAHOBUBINeHCA CKOPOCTbH, 3aMeJICHHA UW TOPMOXKeEHHA. 

Mamepuanoi u memoooi. ba30ii 1a pacueToR HM MOJeMpoOBaHMA CTayia CXeMa J[ByX TpaeKTOpHi MepeMelleHHA V3 
touxw A B TouKy E c yyeToM CHI, 3aTpadeHHbIX Ha 9TH Mporecchl. OntTHMabHoe MepeMellleHve onmpeyemMM c 
TMOMOWbIO upHHuwna IloptasaruHa (TO eCTb ONTHMasIbHoro OpbicTpoyelictBua). I[ponopuMoHasbHoe yiupaBleHve 
IIPHBOJOM TipesCTaBsIeHO KaK MeTO AOCTWKeHHA pesybTaTa. J[y11 WMHHOXOAOBBIX NepeMelleHui MpHBosa TeTabHO 
BH3YasIM3HpOBaHbI (MpeCTaBIeHbl KaK PMCYHKH): CxXeMaTMYecKOe pellieHHe uM pacueTHad cxema. IIpenox%KeHbI 
OPHTHHAIbHbIM CTpyiHbIi WaTaMK C BHYTpeHHeli THEBMaTHYUeCKON CBA3bIO HU THeEBMOMeXaHW4uecKoe JMCKpeTHO- 
TIpOMOpWHOHAasIbHOe YCTpONCcTBO Wa ObicTpoyelicTBHaA KOHTypa yupaByeHua. MatemaTuueckasd MOeuIb BKIOUaeT 
TBYWKCHHe HU TOPMOXKeHHE MOPUWIHA, OamaHc MaCCOBBIX pacxOOB, TaBIeHHe B TOUKaX HM KOHTyp yupaBleHua. Cuctemy 
ypaBHeHu pelwanu Metoqom Pyure — Kytrst B uporpamMMHOM mpogzyKte «Cumuutex» (Simintech). Ilo utoram 


YMCCICTOBAaHuA odo6ujeHHOH MaTemMaTHueckon MOJCIH WOCTpOWJIH 3aBHCHMOCTH HU3MCHCHHA KHHCMaATHUCCKHX, 
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CHJIOBBIX HM IIHe€BMaTHYeCCKHX CBOMCTB TIpHBOJa B peaIbHOM BPeCMeHH IIpH THMOBOM UMKe MO3HIMOHHpOBaHHaA. 
Vndopmauto CyMMUpoBasIN HM IpeCTaBMJIM Kak COBOKYHHOCTb rpadHKos. 

Pe3yismamoi ucciedoeanua. Matematuueckas MOjesIb ctbopMupoBaHa MO KOMIWIeKcy pacueToB. Ona y4HTbIBaeT 
3aBHCHMOCTH, XapaKkTepHble JIA TBWKeHHA NOpWwHA WHeEBMOUMIMHApa. banaHc MaccoBbIX pacxOJOB MccyeqyerTca 10 
ypaBHeHHAM pacxojla ra3a Ip CxKAaTHH B KaMepe, Yepe3 pacipeyeuuTeIM U ApocceNu, B HarHeTaTeIbHOM UM CIMBHOM 
MOJIOCTAX HM B YIIpaBJIAIOWeM ycTpolicTBe. PaccMOTpeHbI HepaBeHCTBA, OMMCbIBaIOLIMe JaBIeHHA B TOUKAX MH KOHTYp 
ylpapienua. ClomwkHad MaTeMaTH4ecKad MOJeJIb pelllaacb B IporpaMMHoii cpeye «CumMuHTex» (Simintech) Metoyqom 
Pyure — KytTpi c u3MeHAeMBIM WaroM MHTerpupoBaHusa. DparMenT paooTsl NporpaMMbl BbIOpaH B KayecTBe OJHOL 
3 WWNOCTpalHi. OH MoKa3bIBaeT, 4TO COT 3afelicTByeT WIA pacueTOB TakKHe [OKa3aTeJIM, KaK: 3ajaHHad 
TIPHBeeHHbIe + KOOPAMHATbI; YHHBepcaylbHadt Ta30Bad MOCTOAHHaA; KOIPMUIMeHTHI %KeECTKOCTH Upy?KHHbI, 
COMpoOTHBIeHHA, ayuaOaTbI WH BA3KOTO TpeHHA B MOpWHe; WaBsleHwe KOMIIpeccopa; Maccy HOABIDKHbIX uacTeli 
MHEBMOMpPHBOa; CHJIy BHCWIHUX COMNPOTHBJICHHH; JHaMeTpbI TpyOonpoBoya, NOPWIHA MHEBMOLWMIMHApa HM TOpMO3HOTO 
yCTpolicTBa; MpOTAKeCHHOCTb XOJa NOPWHA WMIMHApa; MIOWaqWH MOpwHeBbIxX WoMocTeHi WH Apoccenel; WIMHy 
TpyOompoBoya MU ero BHYTpeHHHi oObeM. TakuM oOpa30M, MporpaMMa OlepupyeT 3Ha4HTeJIbHbIM KOMIIJICKCOM 
JAHHBIX, 4TO aeT BOS3MOXKHOCTS MOJYANTb CyU[CCTBCHHbIe HM aj[eKBaTHble pe3ybTaTbl. CxeMaTHYecKH MoKa3aHa 
B3aMMOCBA3b OJIOKOB HM WHarpaMM, HCIIOb30BaHHbIX IPH peleHHH Moyen. Peub ugeT o rpacdbuKax mepememjenui, 
TWIoWaye, TaBweHuli, ckopoctei u Temnepatyp. Ucnomp30BaHbl OOKH C TeEKCTOM IIporpaMMBI HM peqHasHayeHHble 
JUId WHTerpHpoBanna. TakuM oOpa30M MOJy4eHbI MaTeMaTHYeCKad MOJesIb ABTOMATH3MpOBaHHOTO MHeEBMOIpPHBOLa 
TeXHOJOrMYecKorO OOOpyAOBaHHA HU 3ABHCHMOCTH OCHOBHBIX WapaMeTpos ero paootur. [papukKn cBUeTeIBCTBYIOT O 
TOM, 4YTO HCHOMHUTCIBHbIM MCXaHH3M WHEBMONPHBOLa JOJDKHbBIM OOpa30M CylelyeT IpeA102%KCHHOM TpaeKTOpHH. 
O6cyatcdenue u 3akmiouenue. Vitorm padoTbl M03B0JIAIOT PaCCMOTPeTb HECKOJbKO 9TANOB [JIMHHOXOZOBOrO 
TIe¢peMeLIeHHA MIPHBOAa, ONpeAeIUTb BPeEMeHHbIe paMKU 3THX Uporeccos (oT 0 Ao 0,65 cek), a TakxKe (UKCHpyeMBIe B 
JaHHbIC IPOM@XKYTKH H3MCHEHHA aBICHHA HM CKOPOCTH ABWOKeHHA KapeTKH MHeBMOWMIMHApa. Takux 9TaloB MATS: 
pa3roH, JBYWKeHHe C yCTaHOBMBIIelCA CKOPOCTbIO, 3aMeJWICHHe, ABYXKeHHe CO CKOPOCTbIO NO3HUMOHMpoBaHHA 
TopMoxenue. JlaibHeliuMe vccieqoBaHua OyAyT cocpeyoOTOueHbI Ha OMTHMMH3alMH CHCTeMbI JIA coKpallleHHa 


TIPOAOJDKUTCJIBHOCTH HU NOATepKaHiaA TOUHOTO MOSHUMOHHPOBAaHHA TIpH BHCHIHHX BO3J]CHCTBHAX. 


Koroyesbie = CuI0Ba: ~=VIMHHOXOOBOH MHeBMONpHBOT, cCTpyiiHad cCHcTeMa ylpaBJIeHHs, MeXaTpOHHBbIM MOJYJb, 


TIHeEBMaTH4ecKHih AaATANK, MOSHUMOHMpPOBaHHe MHEBMOTIPHBOIa, IpOrpaMMHaad Cpeya «CHMHHTEX>», MeTO], Pyure — KyTrbi 


Baarojapnoctu: aBTOPbI BbIpaxKaroT OlaroqapHOcTE petakWHu UW peleH3eHTaM 3a BHHMAaTeCJIBHOe OTHOMICHHE K CTaTbe 


VW 3aMeudaHHA, KOTOPbIe MO3BOJINJIN MOBbICHTb Ce KadeCTBO. 
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Engineering Research (Rostov-on-Don). 2023;23(3):283—295. https://doi.org/10.23947/2687-1653-2023-23-3-283-295 


Introduction. The performance of the drives is determined by the accuracy of positioning and the speed of 
movement of the coordinates in different operating cycles. Modern fabrication systems are often equipped with 
automatic pneumatic actuators, which are characterized by long-stroke movements. These are, e.g., gantry resistance 
welding machines, coordinate tables, and cutting devices. 

Modern positioning pneumatic actuators for long-stroke movements in fabrication equipment provide a speed of 
up to 30 mm/s and an accuracy of ~1 % of the travel length. Customized actuators provide positioning accuracy of 
0.4 % at speeds up to 100 mm/s. Note that the trajectory of movements is formed by the compressed air flow control 
in pressure or drain pipelines and the pneumatic cylinder sides. In long-stroke drives, the length of such sides reaches 
3 m. Complex thermodynamic processes and compressibility in air flows are the key factors limiting the increase in 
accuracy [1-3]. 

Thus, it is necessary to increase the productivity of the working and engineering processes of the equipment while 


providing the declared accuracy. In this case, it seems promising to create an automated positioning pneumatic actuator 
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for long-stroke movements. The new solution should take into account such characteristics of the pneumatic actuator as 
speed, weight-size parameters, explosionproof, and fire protectability [2, 4]. 

The objective of the work is to obtain a mathematical model and dependences of the key parameters of the proposed 
automated positioning long-stroke pneumatic drive of industrial equipment in the areas of acceleration, steady-speed 
movement, deceleration, and braking. 

Materials and Methods. Figure | shows schematically the transport problem of moving from point A to point E 
along two trajectories. The forces expended on each of the movements are taken into account. The suboptimal 


movement ABCDE (trapezoid) with a simple control algorithm is realized in time 7, >> 7, > min. 


The optimal displacement AFE (bell) was obtained by solving the optimal speed based on the Portnyagin’s principle 


T,, > min, AL < |AL nox | . The result was achieved with more complex proportional control of the drive. The trajectory 


rt 


of movement provided the accuracy of switching motion controls along path L,,. 


Fig. 1. Movement trajectories: 1 — suboptimal ABCDE; 2 — optimal AFE 


In the description of trajectory 1, the switching points are indicated in Latin letters: A — for acceleration of the 


drive; B — for deceleration; C — for positioning speed; D — for stopping. In sections AB and BC, initial acceleration 


and braking are provided up to the positioning speed V,, and further stopping by the braking device AL, < A ae : 
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When moving along the second trajectory, A — switching to acceleration of the drive; F — switching to stop. 

An original jet sensor with an internal pneumatic connection and a pneumo-mechanical discrete-proportional device 
are proposed, which provides increasing the control circuit speed, since feedback in known analogues reduces the 
accuracy of the main engine by about 10-15 % during long strokes [4-6]. 

A schematic solution of a pneumatic positioning drive for long-stroke movements is shown in Figure 2. The drive 
operates in accordance with the suboptimal motion trajectory determined in the optimal speed problem for a given 
positioning accuracy. Here, IIL]1 — rodless pneumatic cylinder of long-stroke movements, which carries out the main 
motion; IIL2, 1113 — brake pneumatic cylinders that fix the drive; CA — jet sensor that determines the coordinate of 
movement, acceleration of the drive, its speed and force; P1 — pneumatic distributor with electropneumatic control, it 
controls the supply to the jet sensor; P2 — main control distributor; P3 — distributor with pneumatic control, it controls 
the operation of pneumatic brake cylinders; [1-4 — mufflers responsible for pressure relief into the atmosphere; 
Jl] — pressure sensor receiving data from the jet sensor; IIJIK — logic controller; LJ] — stepper motor controlling 
the distributor spool; bIIB — air-preparation unit; JJ[P1, J[P2 — throttle with check valve, used to regulate the speed of 


the rodless pneumatic cylinder of long-stroke movements of the main motion [7]. 
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Fig. 2. Circuit design proposal of an automated positioning long-stroke pneumatic actuator 


The drive contains a control system that monitors the position of the carriage of a rodless pneumatic cylinder, slows 
down when approaching the specified coordinates, and sends a signal to the braking gear [1, 4, 7]. 

The speed is increased by the introduction of a discrete proportional control device. Signals generated by the control 
loop are used for control. The device is made in the form of two nozzles and contains compensation measurements. 

Research Results. The design scheme of the pneumatic drive of long-stroke movements is shown in Figure 3. 

Here, IIL{1 — rodless pneumatic cylinder of long-stroke coordinate movements, carrying out the main motion; IL2, 
IIL3 — brake pneumatic cylinders that fix the drive during a stop in the desired position; CA — jet sensor that 
determines the coordinate of movement, acceleration of the drive, its speed and force; P1 — pneumatic distributor with 
electropneumatic control, it controls the supply to the jet sensor; P2 — main control distributor; P3 — distributor with 
pneumatic control, it controls the operation of brake pneumatic cylinders; [1-I'4 — mufflers responsible for pressure 


relief into the atmosphere; J[P1—/[P2 — throttle with a check valve, which serves to regulate the speed of a rodless 
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pneumatic cylinder of long-stroke movements of the main motion; JJ] — pressure sensor receiving data from a jet 
sensor; S — piston area of a rodless pneumatic cylinder of long—stroke coordinate displacements; P1—PS — 
investigated pressures at points 1-5; T1-T5 — investigated temperatures at points 1-5; F:» — friction force in a rodless 
pneumatic cylinder of long-stroke displacements; F,, — viscous friction force in a rodless pneumatic cylinder of long- 
stroke displacements; Fs, — external force in the rodless pneumatic cylinder of long-stroke movements; x — travel of 
the carriage of the rodless pneumatic cylinder of long-stroke movements; V — travel speed of the rodless pneumatic 
cylinder carriage; Crp — spring rate of brake pneumatic cylinders; m — mass to be moved; P; — pressure in brake 
cylinders; Py — pressure in the control channel; f1—/4 — areas of the through sections; Pa — atmospheric pressure; 
d3\1—d33 — diameters of distributor spools; Cappi—Cnpp3 — spring rate of distributors 1-3; xpi—xp4 — displacement of 
distributor spools 1-4; V;1-Vp4 — travel speed of distributor spools 1-4; Fyyi—F's42 — force of the distributor control 


electromagnet 1—2; P,, — compressor pressure; 7, — temperature generated by the compressor. 


Ms p4 


Fig. 3. Diagram of the automated pneumatic drive of long-stroke movements of the industrial equipment 
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The mathematical model was formed with the following assumptions [8-12]: 

— pressure of the compressed air source remained constant over time; 

— thermodynamic process of gas behavior in a pneumatic system was adiabatic; 

— description of pneumatic devices used the ideal gas model, since the pressure of the pneumatic system did not exceed 10 bar; 

— leaks were not taken into account; 

— viscous friction force was proportional to the velocity; 

— expense ratio was recognized experimentally through identification; 

— mass of the moving part was constant; 

— force at the output link of the pneumatic drive was constant. 

1. Equation of motion of the pneumatic cylinder piston [1, 4]: 

=S4(p)-7.)-Fy kay Fug sign a: F, (1) 
Here, S — area of the piston of the discharge and drain sides of the rodless pneumatic cylinder of long-stroke coordinate 
movements of the main motion, m?; P, P, — air pressure in the discharge and drain sides of the pneumatic cylinder, Pa; 


F, 


"ay —— external forces, N; ks: — viscosity friction coefficient; i — friction force, N; o — Boolean parameter: a = 0 at 
Ds < Pon, Ad a= 1 at p, > p,,.,; ps — pressure in the control channel, Pa; p,,— atmospheric pressure, Pa [12-14]; 
m — mass of the moving parts of the drive, kg; Fr — braking force, N. 
Ee a pH . Crp m* (255 ~ Xn ) ? (2) 
where Cnp + — spring rate of the pneumatic cylinder of the brake; ,1 — friction ratio. 
2. Equation of motion of the piston of the braking pneumatic cylinder: 
d?x dx 
m = = Cy m : x m Xin Sin m 1 ee m Kam m : 7 * (3) 
di” (% ) fi dt 
Here, Cup r — spring rate of the brake pneumatic cylinder; xo, — coordinate of the initial compression; S; — effective 
area of the piston of the drain side of the brake pneumatic cylinder, m?; P— air pressure, respectively, in the discharge 


side of the brake pneumatic cylinder, Pa; F's 7— external forces, N; ky: — viscosity friction coefficient. 
3. Mass expenditure balance equations: 


Ga (t)— Gay (1) + Gane (t) = 9: (4) 

G.,,. (t)—G,; (t)-G,, (t) + G.,, (t) =0, (5) 

Gays (t)— Gop (4) — Gyr (1) + Gage (t) = 0, (6) 

Giys (t)— Gaps (4) — Gaya (t) + Gag (1) = 0, (7) 

Cry (t)— Gps (4) + Ge (1) = 0 - (8) 

Here, G,,(¢) — gas mass flow rate under compression in the chamber; G,,(t),G,,(¢)— mass flow through 


distributors; G 


mi (t) and G,,,.(¢) mass flow rate in the discharge and drain chambers of a rodless pneumatic cylinder of 


long-stroke movements; G,,(t), G,,,(t), G,,.(t), G,,3(t), G,,(t) — mass flow through throttles in the drain line, at the 
inlet to the nozzle unit of the jet sensor and at the outlet of the nozzle unit; G,,,(¢), G,,.(4), G,, (0, G,,.4(4) — mass flow 


rate in the control channels of the control device, at the outlet of the control device, distributor of brake pneumatic 
cylinders [7]. 
awViod 
(=P ) 
4-E-R-T dt’ 
where p= 1.288 kg/m? — air density; V — volume of the chamber; R = 287 J/(kg-k) — gas constant; E — volume 


modulus of elasticity of air; T — temperature at the point; “ — pressure variation at the point [6]. 
t 
W, dp 
G,,, (t)=— - + 10 
mi (1) = k-R-T dt’ 
W, dp 
G,,,. (t) =—+—: =, 11 
mo (t) k-R-T dt oo 
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where W, and W, — current volumes in the pressure and drain sides of the base pneumatic cylinder, m+; & — adiabatic 
exponent (for air, k= 1.4). 
4. Pressure equations in points: 


ap, BART [p,, oe {2 va Top? PL a 12 
dt S(x+Xq, )JE, (x+X, ae ee ae 


k-1 


k- f.-.JR-T. 2k F 
Ps fa BT = Se ae (13) 


ie. dt” 


dt S(L-x+Xy) JE, P, 

dp, hv EST (2) ra 
= W “VP TE a ae rae (14) 
dp, _ k- | R T, k- Ss: R T, Ps —— 

PE MEE 


kK: R-T, 
dps _ ff = aA 


dt W,- 


Here, k — adiabatic exponent; R — gas constant, J/ — ; Tw — air temperature in the main line, K; Pa — 


atmospheric pressure, Pa; P; — Ps — pressure in the flow parts of the pipeline, Pa; W; — Ws — volumes of the flow 
parts, m?; €;— &7 — resistance values in the line; fi —_fs— areas of the flow sections of the pipeline, m’; L — maximum 
stroke of the piston, m; xo1, x02 — ratio of the initial volumes of the pneumatic actuator to the useful area of the piston 


side of the pneumatic cylinder, m; — piston speed, m/s. 
t 


5. Control loop equations. 
Equation of motion of the distributor spool 1: 


M, » i Cpa ta Fa sign( Fi (2 Fa (17) 

Equation of motion of the distributor spool 2: 
pp = eg Sn —P -sign( Se Fs (= Fyz: (18) 

Equation of motion of the distributor spool 3: 
ea 

Equation of motion of the distributor spool 4: 


Here, S,,— area of the distributor spool end face, m*; P3, Ps, Py — pressure in the control channels, Pa; F’.— resistance 
forces to the movement of the distributor spool, N; Fy — reaction forces of stops, N; F'sm — force of the electromagnet 
acting on the distributor spool, N; cup p — spring compression ratio, N/m; m;, — weight of the distributor spool, kg. 

To solve the mathematical model, SimInTech software product was used. We applied the Runge—Kutta method with 


a variable integration step (Fig. 4—5). 
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Fig. 4. Model fragment in SimInTech program 


Here, P,1—P,3 — blocks with the program text; f1, f2 — graphs of the output of areas | and 2 of the jet sensor 
throttle; fis: — block of the problem of the areas of the flow sections of the jet sensor; fis2 — graphs of the output of the 
areas of the flow sections of the jet sensor; d71-d7T8 — integration blocks; P1—P4 — graphs of the output of the 
received pressures at points 1-4; P — general schedule for the output of all pressures; V — graph of the obtained speed 
output of the rodless pneumatic cylinder carriage; x — graph of the received movement output of the rodless pneumatic 
cylinder carriage; x, V — general graph of the output of the received movement and the speed of the rodless pneumatic 
cylinder carriage; 7 — graphs of the output of the obtained temperatures; ¥3 — graph of the received movement output 
of the control distributor spool. 

2 Se Re nperpagoo an L=0.43 
ee S=(pi*D*2)/4; 


QA *eOx|1|Se/At|E FUR elrv : 
1 input pl, p2, p3, p4, V, x,fs1, fs2, ds; St=(pi*Dt*2)/45 


output pl, p2, p3, p4, V, x, d, T1, 12, T3, T4, TS, T6,p5; | kt=0.4 
init p1=100000,p2=100000, p3=100000, p4=100000, V=0, x=0; begin 
xk=0.31 30 | if x<xk+@.0577 then 

- R=287; 16=0 
Tm=293; foe 
ks1=653 else 
ks2=100; p6=pm 
ks3=70; end 

10 ks4=70; 5 begin 
d3=0.001; ; . 
d5-0.006; if x>=0.363 then 
pm=6*1045; p5=pm 
pa=1.01325*10°5; else 
k=1.45 p5=0 
x01=0.008 ae | end 
02=0.008; 

Acs N=p5*St; 
=6 5 
Fvn=20; Ft=kt*N; 

20 dtr1=0.005; begin 
dtr2=0.005; if x<xk then 
peel d2=0.008 
Dt=0.04; 1 
kvt=320; bag 

d2=0.0004 
end; 


b) 
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“begin a 
: begin 
so | if x<xk+@.0577 then 4f p3>=600000 then p3=600000; 
d=ds = 
1 5 begin 
i 7 a lls then p3=pa; 
— = en 
s 20 begin 
end; if p4<=pa then p4=pa; 
© begin end 
if x<xk then T= (p/p) *( (K-27) *Tm5 
5-0 T2=(p2/pm)*((k-1)/k)*Tm5 
a T3=(pm/pm)*((k-1)/k)*Tm; 
else La ila a 
i T5=(pm/pm)*((k-1)/k)*Tm; 
5=(pi*d5*2)/4 T6=(p4/pm)*((k-1) /k)*Tm5 
60 end; G11=(k*f1*sqrt(R*Tm) /(S*(x+x01)*sqrt(ks1) ))*sqrt(abs(pm*2-p1*2) ); 
1=(pi*d*2)/4; 20 G12=(k*f3*sqrt (R*T3) /(S*(x+x@1)*sqrt(ks2)))*(p3/p1)*((k-1)/(2*k))*sqrt(abs(p1*2-p3*2)); 
- - G13=(k*f3*sqrt(R*T5) /(S*(x+x01)*sqrt(ks2) ))*(p4/p1)*((k-1)/(2*k))*sqrt(abs(p1*2-p4*2)); 
#2=(pi*d2*2)/4; G14=(k*p1/(x+x01))*V; 
£3=(pi*d3*2)/4; G21=(k*f2*sqrt (R*T2) /(S*(L-x+x02)*sqrt(ks2) ))*(p2/pa)*((k-1)/(2*k))*sqrt(abs(p2*2-pa%2)); 
G22=(k*p2/(L-x+x@2) )*V; 
Ltr=0.01; G31=(k*#3*sqrt (R*T3) /(W1*sqrt(ks3)))*sqrt (p1*2-p342); 
cpr=500000; G32=(k*fs1*sqrt(R*T4) /(W2*sqrt(ks4)))*((p3/pa)*((k-1)/(2*k)))*sqrt(abs(p3*2-pa*2)); 
ee eee G41=(k*f3*sqrt (R*T5) /(W1*sqrt(ks3)))*sqrt(p1*2-p4%2) ; 
W2=Ltr*((pi*dtr2*2)/4); G42=(k*fs2*sqrt(R*T6) /(W2*sqrt(ks4)))*((p4/pa)*((k-1)/(2*k)))*sqrt(abs(p4*2-pa*2)) 5 
iS begin p1'=G11-G12-G13-G14; 
ey 100 p2"=-G21+G622; 
if L<x then Fy=(cpr*(x-L)) p3"=G31-G32; 
70 | else p4'=G41-G42; 
if (L>=x) or (x>=0) then Fy=0 ae ‘io ogeticy haciad epee 
else if x<@ then Fy=cpr*x SR a ea rr mee, fe nea eye rahe 
c) d) 


Fig. 5. Programming block in SimInTech: a — part of the source data; b — part of the logical functions; 


c — part of the assignment of variables; d — skeletal code 


The study of the generalized mathematical model of the proposed drive made it possible to obtain graphs of the 
behavior of the drive during acceleration, deceleration, and stop (Fig. 6), describing the changes in kinematic, power 


and pneumatic properties of the drive during a typical positioning cycle in real time [15]. 


P5, Pa V,m\s x,m P4,Pa P3,Pa P2,Pa P1,Pa V,m/s P3,Pa P4,Pa x,m P2,Pa Pl, Pa P5, Pa 


600000 | 9.46 1 460000 | 480000 


580000 580000 1 0000 580000 
560000 | 0.34] 9-95 | 560000 460000 
540000 0.9 ]540000 {420000 
s20000 | 82 520000 { 4o0000 {42° }640000 
500000 0.3| 851500000 | senqqq {420000 4620000 
480000 0.8 |480000 500000 
460000 | 0.28 460000 {960000 480000 
0.75 
_— 026 pas 340000 | 380000 |4go900 
0000 0.7 {420000 | 320000 440000 
400000 | 0,24 400000 os 
0.65 | sa999 4300000 1420000 i 


380000 340000 Aaa : 
360000 | 0.22 0.6 |360000 |280000 1400000 i | lil HH : 
940000 | 45] .55, {940000 | 260000 320000 1550000 | i 
320000 “1320000 1.9000 1300000 }360000 i | H 
300000 0.18 0.5 {300000 140000 i 
280000 280000 |220000 {280000 7 
260000 | 9:16} 45 1260000 4200000 | agqqqq 42200 i 
' 

240000} 9 44 o4 — 180000 | soon {300000 : 
220000 0.35 20000 1 1280000 ‘ 
200000} 0.12 200000 | 160000 1260000 3 
220000 7 

180000 0.3 | 180000 | 140000 5 i 
160000 01 160000 | 120000 | 200000 {24200 i 
140000 0.25 | 440000 220000 i 
1 

i 

' 

i 

i 

; 


0.08 00000 
420000 0.2] 120000 |" 180000 t>o9000 val 
100000 0.06] , ,,. {100000 80000 | 1g9000 1180000 
80000 “"°1 go000 | 60000 pee 
60000 | 9.04 0.1} 60000 140000 
40000 140000 

40000} go 5] 49000 120000 
20000 0 20000 | 20000 120000 

o 0 Y) 0 | 190000 100000 


0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 045 0.5 0.55 0.6 0.65 0.7 t,sec 
Fig. 6. A set of graphs based on the data of a generalized mathematical model 


Discussion and Conclusion. The graph shows a long-stroke movement along the trajectory proposed in Figure 1. 
The operation of the drive consists of several stages. 

1. 0-0.1 sec. Acceleration. The pressure in the pressure side increases to 4 bar, the pneumatic cylinder carriage 
speed — | m/s. 

2. 0.10.38 sec. Movement at a steady speed. Pressure in the discharge side — about 3.8 bar. Pressure in the drain 


line increases to 1.6 bar, the pneumatic cylinder carriage speed — 0.85 m/s. 


http://vestnik-donstu.ru 


3. 0.38-0.5 sec. Slowing down. Pressure in the pressure and drain sides increases. The pneumatic cylinder carriage 
speed decreases to 0.075 m/s. 
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4. 0.5-0.65 sec. Movement with the speed of positioning. Pressure in the pressure and drain sides increases, the 
pneumatic cylinder carriage speed — 0.075 m/s. 

5. 0.65 sec. Switching to braking, activation of an external braking gear. 

The graphs obtained confirm that the long-stroke movements of the pneumatic actuator are performed in accordance 
with the proposed trajectory (Fig. 1), and the control system functions properly. Further research will focus on 


optimizing the system to reduce the duration and maintain accurate positioning under external influences. 
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3aABNeHHbIU BKIAO Coaemopos: 

JJ.A. Kopotsra — moyroropka TekcTa, cocTaBileHHe MaTeMaTHYecKOl MOJeIM, MOJeMpOBaHHe B IIporpamMe, 
aHalIH3 pe3yJIbTaTOB HccIeqOBaHHA. 

B.C. CunopeHko — HaydHoe pyKoOBOACTBO, POPMHPOBAaHHe OCHOBHOM KOHIICHIMH, IesIM MW 3ayad MCCIeqOBaHHA. 


CI. WpuxoypKo — ocdbopmienne 4 JopadoTKa TeKCTa, KOPpeCKTHPOBKa BbIBOJOB. 
Kongauxm unmepecoe: aBropbl 3atBIAIOT OO OTCYTCTBHH KOH@JIMKTa HHTepecos. 


Bce aemopol npowumaau u odobpusu OKOHYaMeNbHoIU BapuaHm pyKoNucu. 
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